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Summary 

The practical methods of measuring reverberation time are discussed briefly. 
The design of a microprocessor controlled system for measuring reverberation time 
in a wide range of different types of rooms is described. The system employs 
analogue signal processing to provide the loudspeaker excitation signal and to 
receive and process the microphone output signal. The sequences of operations 
are controlled by a microprocessor which also controls the display of the sound 
pressure level as a function of time using a conventional television format monitor. 
The microprocessor is also programmed to calculate and display the numerical 
value of the reverberation time. Measurements can be carried out in any of 24, 
1/3-octave wide frequency bands from 50 Hz to 10 kHz. An interface to an agreed 
international standard is provided to enable remote control of the instrument 
as part of a larger measurement system. Field trials, using the prototype instrument, 
have shown that it substantially fulfills all of the design requirements with only 
minor faults and modifications. 
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1. 



Introduction 



The control of the acoustic environment is 
an important part of the design of many rooms. 
This requirement is particularly evident in a 
broadcasting organisation which uses a large 
number of studios, control rooms, listening rooms 
and film dubbing theatres, all of which require that 
the internal acoustic environment is 
well-controlled. One simple measure of this, which 
has been in use for many years, is the 
"reverberation-time,", which is defined as the time 
taken for an initial sound pressure level to decay 
by 60 dB after the removal of the sound source. 
In general, the rate of decay of the sound energy 
will be a function of frequency and, as such, 
should be specified and measured at a number of 
different frequencies. For broadcasting purposes, 
it is usual to design each room so that the 
reverberation time is the same at all frequencies 
within the range normally considered, usually 
50 Hz to 8 kHz or 10 kHz. Nevertheless, it is 
still necessary to carry out the measurements at 
each frequency to confirm that the design 
specification has been achieved. Although a 
number of different methods have been proposed 
for measuring reverberation times 1_4 , the original 
method and the one still most commonly used 
consists essentially of exciting the room with a 
sound source and recording the time function of 
the decay of the sound pressure level after the 
sound source is removed. From the recording, the 
rate of decay of the sound pressure level in decibels 
per second can be measured and the reverberation 
time calculated. For each measurement frequency, 
this procedure is repeated for a number of decays 
to reduce the effects of random variations, and at a 
number of microphone positions to obtain the 
average reverberation time throughout the volume 
of the room. In an average size of room this 
measurement would be carried out at five different 
microphone positions and at 24 different 
frequencies spaced at 1/3 -octave intervals. An 
average of approximately 3 decays would be 
recorded for each position and frequency. Thus a 
total of approximately 360 decays would be 
recorded and analysed. 

This measurement, if carried out manually 
using a pen-recorder to record the decays, is 
extremely time-consuming and laborious, 
especially as, in the BBC, it must be repeated 



for about 50 — 60 new or refurbished rooms per 
year. Therefore, earlier attempts had been made 5 - 6,7 
to give the operator a direct-reading of the 
reverberation time. The most successful of these 8 ' 9 
displayed the logarithm of the decay of the sound 
pressure level on a cathode-ray tube with a 
long-persistence phosphor. An adjustable graticule, 
which was calibrated directly in reverberation time 
and which could be manually aligned with the 
average slope of the displayed decay, was provided 
to give a direct reading of reverberation time. 
Using this equipment, the average of several decays 
could be estimated visually and the corresponding 
reverberation time obtained from the graticule. 
The settings of the sound-source frequency and 
level, the necessary filter centre frequency and 
the microphone amplifier gain were still carried out 
manually by the operator. Despite this, the 
equipment greatly simplified the measurement of 
reverberation time and, with later improvements, 
has been in constant use for about 30 years. 
Although it is still serviceable, some of the spare 
parts arc becoming difficult to obtain. Also, 
by modern standards the complete set of equip- 
ment including the sound source, filter set, micro- 
phone amplifier, and calibrator is very bulky and 
heavy. 

In 1969, an alternative method was 
developed which was based on the analysis of tape 
recordings by a fixed, large mainframe computer. 
A tape recording with a prerecorded excitation 
signal was replayed and the response of the room 
to this signal was recorded on a second tape. This 
second tape was subsequently returned and 
analysed by the computer. Although this method 
did not give results immediately, it was used for a 
number of years and was particularly suitable for 
obtaining measurements using non-specialist 
operators. The excitation tape could be sent by 
post and the measurements carried out by the 
local staff. 



The advent of modern circuit elements has 
since made it possible to carry' out the same 
functions as the original portable equipment at 
a very much lower cost and weight. In addition, 
the measurement procedure is well defined and 
is eminently suitable for automation by the 
incorporation of a microprocessor-based control 
system. 
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2. 



Choice of the measurement method 



A number of possible methods for measuring 
reverberation time have been proposed 1 " 11 and 
some commercial measuring systems are available. 
These can be divided into two major categories, 
the first based on the use of an impulse of sound to 
excite the room. Subsequent computer processing 
of the microphone signal extracts the reverberation 
time as a function of frequency. These methods 
generally suffer from the disadvantages of poor 
signal to noise ratio unless many pulses are 
averaged and the requirement for a comparatively 
powerful computing system. The second major 
category consists of those methods which are 
essentially the same as the basic method described 
above in the introduction, followed by very much 
less complicated data processing to extract the 
reverberation time. This second category has one 
major advantage in that, at least in principle, the 
shape of the sound pressure level decay as a 
function of time can be displayed directly as it 
occurs. This facility is valuable in the diagnosis 
of certain types of faults in the acoustic 
environment. A further major advantage of this 
second category is that the component parts of 
the system, such as the signal source and the 
frequency band selection filters, are identifiable 
separate elements which could be applied in other 
types of acoustic measurements. 

Three main considerations led to the 
decision to design the measuring system rather 
than to buy commercial equipment. The first of 
these was that no commercially available 
equipment satisfied the requirements completely, 
although it would have been possible to use any of 
them by accepting some restrictions and by 
modifying the existing measurement procedure.* 
The second consideration was that of total cost. 
It would have been a serious limitation to the 
normal programme of work to purchase only one 
set of equipment. A cost analysis showed that, 
provided that the development costs could be 
shared between at least two sets of equipment, it 
was more economical to design and build the 
equipment than to buy it. Additional benefits 
obtained would be the supply of equipment to 
other departments within the BBC at a very low 
unit cost and the possibility of recovering some of 
the design costs by licensing agreements with 
interested commercial manufacturers. The third 
consideration was the possibility of using the 

* One main requirement was for the display of the time function of 
sound pressure level. Some commercial equipment had this facility 
but was very slow and required an additional display device 
(oscilloscope). The remainder did not have the facility. 



design of the component elements for other 
purposes, something which would be difficult (and 
unethical) with commercial equipment. (See 
Section 9). 

3. System outline and operating principles** 

Fig. 1 shows a block diagram of the system. 
It consists of a source of the signal for exciting the 
room, a microphone amplifier, a filter to reduce 
the effects of noise at frequencies other than the 
measurement frequency, and a unit to convert 
the amplified microphone output to a signal 
representing sound pressure level. These units 
are controlled by a microcomputer which generates 
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Fig. 1 - System block diagram. 

**The instrument described in this report is similar to that described 
in reference 11, the description of which was published after the 
work described in this report had been completed. 
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the appropriate control and timing signals. A 
display unit to show the decaying sound pressure 
level as a function of time and a keyboard to 
control the system operation are also included. 

For the storage and subsequent processing 
of the sound pressure level function, the signal 
must be converted from analogue to digital form at 
some stage. A rough cost-benefit analysis for 
different possible locations of the analogue to 
digital convertor is given in an Appendix (Section 
12.1). It was decided that the conversion from 
analogue to digital form should be done after the 
conversion of the microphone signal to the form 
of the logarithm of the mean sound-pressure 
level. This meant that all of the audio-frequency 
circuits were analogue and were very much easier 
to implement for this application than their 
equivalent digital counterparts. 

In operation, the microprocessor interprets 
inputs from the keyboard and sets the operating 
frequency and reverberation time range. 
Reverberation time measurement commences on 
receipt of a 'start' instruction. Each measurement 
cycle is divided into two parts, an excitation period 
and a decay period. When the start instruction 
is received, the output of the excitation generator 
is enabled for a time period which is long enough 
to allow the reverberant sound energy in the room 
to reach its asymptotic value. This time period 
depends on the reverberation-time range to which 
the instrument has been set. During this period, 
the sound pressure level received by the 
microphone is filtered and sampled to enable the 
microphone amplifier gain to be set for the 
subsequent decay period. Near to the end of the 
excitation period the microphone amplifier gain is 
set to a value calculated to raise the microphone 
signal to a suitable level for subsequent processing. 
The short delay between the setting of the gain 
and the termination of the excitation output 
allows the microphone amplifier circuits to settle. 

At the end of the excitation period, the 
output is disabled and, simultaneously, the display 
of the sampled sound pressure level as a function 
of time is started. In addition to being displayed, 
the sound pressure level data is also stored for later 
computation. This cycle of excitation and decay 
can be repeated automatically. 

The system as described is essentially a 
modern replacement of the equipment which has 
been in use in the BBC for many years. However, 
the incorporation of the extended capabilities of a 
microprocessor permitted the addition of further 
features, not present in the old equipment, at no 



additional component cost. Most important of 
these is the ability to generate the cumulative 
average of the data from several decays (at the 
same frequency), thus reducing the effect of the 
large random variations of the sound pressure level 
time function common to this type of 
measurement, particularly at low frequencies 12 . 
This average time-function is also displayed. In 
addition, from the stored data representing the 
average decay function, the average decay rate in 
dB/s and, hence the reverberation time is calcu- 
lated arithmetically, as described in Section 7. 

4. Design specifications 

4.1. Frequency range 

For normal room acoustic applications, a 
measurement frequency range is required which in- 
cludes most of the audible frequencies. The freq- 
uency range is also required to be sub-divided into 
parts small enough to give reasonable resolution of 
the reverberation time as a function of frequency. It 
was decided that a range of 50 Hz to 10 kHz, 
subdivided into 24, 1/3 -octave wide bands would 
be adequate to cover most requirements.* 

4.2. Reverberation time range 

The reverberation times encountered in 
broadcasting studios and other areas range from 
about 0.15 seconds in small, heavily-treated rooms 
to about 15 seconds in large orchestral studios and 
concert halls prior to the installation of the 
acoustic treatment. The measuring instrument has 
to be capable of giving accurate results over this 
range of reverberation times. To achieve this, 10 
ranges were thought to be desirable. The shortest 
range, to 0.5 seconds, was limited by the ability 
of the microprocessor to process the data. This 
range would give meaningful measurements of 
reverberation time down to about 0.12 seconds. 
The longest range, to 15 seconds, would be 
capable of giving accurate measurements up to a 
reverberation time of 15 seconds and, by utilising 
less than the full dynamic range, somewhat less 
accurate results up to at least 30 seconds. There 
was, however, no limitation to the longest range 
which could have been incorporated. A time 
resolution of 200 samples for the full scale range 
was set by the display limitations and this was 
thought to be adequate. 

4.3. Signal amplitude range 

The maximum usable amplitude range of 

* Routine measurements in the BBC had previously been made only 
in 16 half— octave wide bands over the same frequency range. 
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microphone signals in this type of measurement is 
limited at one extreme by the maximum available 
sound power output of the loudspeaker and at the 
other by the background noise level in the room 
being measured. In practice, although dynamic 
ranges of 80 dB can be achieved in the 
mid-frequency range, it is rarely necessary to use 
such a large range. A range of about 60 dB with a 
resolution of %dB (240 discrete levels) was 
considered adequate. This fits neatly and 
economically with commonly available, 8 bit 
microprocessors which give 256 discrete levels for 
one byte. This dynamic range limitation is one of 
display and computation only — the absolute range 
of input signals is set by the gain -controlled 
microphone amplifier as described in Section 3. 

5. Detailed description of the individual 

units ' ' ' : ' ' . 

5.1. General •■■■■•■■ _ 

This section (and Section 6) inevitably con- 
tain terms and expressions common in the descrip- 
tion of electronic and ' microcomputer equipment 
but less so in acoustics. A short glossary of these 
terms is given in an Appendix (Section 12.2). 

5.2. Excitation generator 

The excitation generator unit is required to 
provide the signal for driving the loudspeaker. The 
signal source was chosen to be band-pass filtered 
random (actually pseudo-random) noise as, of all 
the suitable types of signal, this one is compara- 
tively easy to generate. The output signal 
bandwidth defines the measurement frequency and 
has therefore to be reasonably accurately 
controlled. Fig. 2 shows a simplified circuit 
diagram of the excitation generator. The design of 
the unit is based on a single-chip, integrated 
bandpass filter (Reticon R5604). This integrated 
circuit provides, by means of a network of 
switched capacitors and amplifiers", a filter with a 
six-pole Chebyshev response. The pass-band centre 
frequency is proportional to the clock frequency 
and the constant fractional bandwidth is equal to 
1/3 -octave. The input signal to the filter is derived 
from a 24-stage, maximal-length feedback shift 
register. The clock frequency for the filter and 
shift register is provided by a crystal oscillator 
via a chain of seven binary division stages. By 
selecting succesive outputs of the divider chain, 
the centre frequency can be set at octave intervals 
in the range 8 kHz to 63 Hz. Unfortunately, 
there is no elementary, accurate relationship 
which will subdivide an octave into 1/3-octave 
intervals. The simplest solution was to provide 



two additional crystal oscillators with frequencies 
spaced by ± 1/3-octave relative to the centre 
frequency oscillator. By selecting the' appropriate 
oscillator as the input to binary divider chain, 
the required 1/3-octave centre frequency is realised. 
The complete output frequency range is thus 50 Hz 
to 10 kHz in 1/3-octave steps. 

The output signal from the filter contains 
a low-level residue of the clock frequency and, 
because of the sampled nature of the signal 
within the filter, spurious outputs at multiples 
of the clock frequency ± centre 'frequency;' As'the 
filter clock 1 frequency is equal to 54 'times-" -the 
centre frequency, these spurious outputs can be 
reduced to an insignificant (but hot inaudible 1 ) 
level by a relatively simple, two-pole analogue 
low-pass filter: This simple filter was designed to 
be switchable only in octave increments ' and the 
cut-off frequency was set to be approximately 
one octave above the output frequency. The 
rejection of the spurious outputs is thus greater 
than 50 dB." : - " ' "' " ■ '" " - 

The output level from the excitation unit is 
controlled by ' ah integrated, digitally-controlled 
attenuator (Analogue Devices, type AD7110). The 
attenuation can be controlled in 59 steps of 1 : .5 dB 
from OdB to 88.5 dB, plus an additional step to 
give infinite attenuation. A provision has also been 
made for using an external signal source instead of 
the internally generated, band-pass filtered random 
noise. 

To control the unit, an interface to the 
microprocessor address and data buses allows the 
unit to be used as a peripheral output device on the 
microprocessor system. The operating frequency 
and output level can be set by writing, from the 
microprocessor, a single 8 bit binary word into the 
two addresses controlling the centre frequency and 
output level respectively. • 

5.3. Microphone amplifier and filter 

The microphone amplifier was designed to 
accept the output from normal moving-coil 
microphones at an impedance level of -approxi- 
mately 200 — 300 ohms. It was specifically designed 
to accept the output from an STC type 4021 H 
microphone which has been proved ' to be 
sufficiently sensitive and robust for this application 
but it will work with any electromagnetic 
microphone which has a sensitivity of about— 53 dB 
relative to lV/Pa including the microphone 
transformer. " No provision was made in the 
prototype for using capacitor microphones directly 
or for providing phantom power for any type of 
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Fig. 2 - Simplified circuit diagram of excitation generator. 



microphone. To eliminate the requirement for 
manual setting of the microphone gain, the 
amplifier was designed to be gain-controlled by the 
microprocessor system over a range of gain from 
40 dB to 100 dB. This is sufficient range to cover 
the majority of applications. With the type 402 1H 



microphone, the lowest gain setting (+40 dB) is 
such that the amplifier will overload at sound 
pressure levels of about 110 — 115 dB re 20juPa. 
With the highest gain setting (+100 dB), the range 
of sound pressure levels displayed by the complete 
instrument is sufficient to include the theoretical 
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intrinsic noise voltage of a 300 ohm source 
impedance. To maximise the usable range, the 
amplifier was designed to have as low a noise 
figure as was possible without unduly compli- 
cating the design. 

To achieve the design requirements, the 
amplifier was designed with two stages. The 
first or preamplifier stage consisted of a 
conventional, low-noise transistor amplifier with an 
input impedance-matching transformer and a fixed 
voltage gain of +40 dB. The typical noise figure 
for this stage was 1 dB for a 300 ohm source 
impedance. The second or main stage of the 
microphone amplifier was based on the same 
type of digitally controlled attenuator as was used 
in the excitation generator, but arranged within the 
feedback loop of an integrated circuit amplifier 
to give gain instead of attenuation. A total range 
of gain from to +60 dB is available from the 
second stage. The comparatively poor noise 
performance of the integrated circuit amplifier is 
not significant because of the preceding 40 dB 
gain in the first stage. 



A band-pass filter was included following 
the microphone amplifier to improve the signal-to- 
noise ratio in the frequency range being measured. 
One requirement was that the filter should have a 
fast response in the time domain to allow short 
reverberation times to be measured by the 
complete instrument. The time-domain response 
of the integrated filter used in the excitation 
generator was unacceptable. In addition, because 
the time response of a linear-phase filter is 
inversely proportional to the bandwidth, a 
1/3-octave filter bandwidth was unacceptable for 
the lowest frequency octave (50, 63 and 80 Hz) if 
the requirement to measure reverberation times 
down to 0.15 seconds was to be met. A simple (in 
principle) four-pole, band-pass filter with a one 
octave bandwidth was acceptable. The basis of the 
filter design is shown in Fig. 3 . By switching each 
of the four variable resistors to one of 24 different 
values, the centre frequency of the filter could be 
set in increments of 1/3-octave steps. Such a 
switchable filter occupies a large area of printed 
circuit board and uses a large number of electronic 
switches and resistors. It is, however, still 
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comparatively cheap and reliable. Unlike some 
other types of active band-pass filters, neither the 
gain nor the bandwidth is seriously affected by 
component value tolerances. It can be 

satisfactorily constructed using readily available 
2% tolerance resistors. 



An interface to the microprocessor 
address and data signals allows the microphone 
amplifier gain and filter centre frequency to be 
controlled in the same way as the excitation 
generator. 
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of-53 dBre lV/Pa). 



A fundamental problem was not 
encountered until the completed prototype had 
been assembled. The proximity of the internal 
television picture monitor and the high-gain 
microphone amplifier resulted in line-scan 
frequency components being induced into the 
low-level microphone circuits. No reasonable 
amount of shielding was sufficient totally to 
overcome this problem and so a low-pass filter 
was designed to reject these line-scan frequency 
components. This filter has a pass-band with 
0.5 dB ripple extending to 11.3 kHz and a 
minimum stop-band attenuation of 47 dB. An 
elliptic filter was used and the design is such that 
the first null (>60 dB attenuation) in the filter 
response is at 15.625 kHz. (line-scan frequency). 
The harmonics of the line-scan frequency are 
much less significant because of the poor 
frequency response of the succeeding stages above 
about 30 — 40 kHz. Also, the filter attenuation of 
these harmonics is not less than 47 dB. The filter 
was placed at the output of the microphone 
amplifier board after the switchable band-pass 
filter. 



Some in-band pickup remained from the 
harmonics of the monitor field-scan frequency 
(50 Hz) and this limited the noise performance 
of the microphone amplifier at frequencies below 
about 300 Hz. This is not a significant limitation 
because most rooms have moderately high levels 
of low-frequency noise, especially if mechanically 
ventilated. Fig. 4 shows the equivalent sound 
pressure level of the inherent noise level of the 
completed prototype instrument for a microphone 
sensitivity of -53 dB re lV/Pa (as for STC 
4021H) and, for comparison, the theoretical 
equivalent noise level of a 300 £2 resistor at 293 K 
in a one octave-wide measurement band. In 
practice, the microphone amplifier filter has an 
effective noise bandwidth somewhat greater than 
one octave. 



Measured system noise level. 

Theoretical equivalent noise 
level of 300 £2 resistor at 
293 k in octave bandwidth. 



5.4. Rectifier, integrator and analogue-to- 
digital convertor 

The microphone amplifier and filter output 
signal has to be converted to a signal representing 
the sound pressure level. This requires the signal 
to be rectified and then filtered to remove the 
rectified signal frequency components. A precision 
rectifier followed by a single pole, low-pass filter 
(integrator) is sufficient. Unlike in the earlier 
equipment referred to in Section 1, the rejection 
of the rectified signal frequency components is 
not required to be so high because of the electronic 
trace averaging facility included in the new 
equipment. Any residual signal frequency 
components in the individual traces tend to be 
reduced during successive traces by the averaging 
process. This allows the integrator response time 
to be somewhat shorter than in the earlier 
equipment 8 ' 9 , thus allowing shorter reverberation 
times to be measured. The final design has an 
integrator response time equivalent to a 
reverberation time of approximately 0.065 seconds. 



The output from the low-pass filter 
represents the short term average of the magnitude 
of the band-pass filtered microphone signal, 
whereas the sound pressure level is defined as the 
root mean square of the filtered microphone signal. 
However, the signal present at the output of the 
band-pass filter is narrow-band, filtered random 
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Fig. 5 - Simplified circuit diagram of rectifier and logarithmic law analogue-to-digital convenor. 



noise, further filtered by the time response of the 
room being measured. Except for a constant 
scaling factor, the arithmetic average of the 
magnitude of this signal is the same as the root 
mean square to a sufficient degree of accuracy for 
the purpose of measuring reverberation time. 
A worst-case calculation based on a single, centre 
frequency sinusoidal waveform, amplitude 
modulated to a depth of 50% by a lower frequency 
sinusoidal waveform gave an error of less than 
0.5 dB in taking the scaled mean instead of the 
true r.m.s. In practice, with the effective 
modulation consisting of band-limited random 
noise to a lower modulation depth than 50%, the 
errors will be much less than 0.5 dB. 

The output from the low-pass filter is 
converted to a digital form by simultaneously 
starting a ramp waveform and a binary-coded 
digital counter. The counter is stopped when the 
ramp waveform becomes equal to the input signal. 
In this way, a digital representation of the 



amplitude of the input signal is generated in the 
counter. By making the ramp waveform 

exponential rather than the more usual linear 
ramp, a direct conversion to a logarithmic (i.e. 
decibel scale) is obtained. 

Over a working range equal to 64 dB 
(= 1600 : 1 voltage range), the inherent input 
voltage offsets of the rectifier and the comparator 
amplifiers would have presented serious difficulties. 
To overcome these, the entire system is duplicated 
(See Fig. 5) with one half being preceded by an 
amplifier with an a.c. gain of 32 dB. Input signals 
in the range to 31.75 dB are converted directly 
to a 7-bit digital form as described above. Input 
signals in the range —32 dB to —63.75 dB are 
amplified by 32 dB and then converted. Whichever 
of the two conversions proceeds to completion 
controls the most significant bit of the 8-bit digital 
output representing the logarithm of the amplitude 
of the input. In this way, the effects of the input 
offset voltages of the active components are 
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are reduced by the equivalent of 32 dB. The 
resulting complete unit requires only two adjustable 
components, to set the gain of the 32 dB amplifier 
and the time constant of the exponential ramp 
waveform. This method is more economical 
and stable than using selected low offset devices or 
offset-trimming components. The resolution of 
the digital output is 64 dB / 2 8 , that is % dB per 
step. 

An interface to the microprocessor address 
and data signals is provided to enable the transfer 
of the data to the microprocessor. The data 
sampling and conversion, which takes a maximum 
of about 1280 jusec, is initiated by the reading of 
the previous data. 

5.5. Display 

The display is required to show both 
alphanumeric data and the graphical data 
representing the individual and the averaged decay 
of the sound pressure level as functions of time. 
A resolution of 256 x 256 picture elements allows 
the two graphical outputs each to be displayed 
with an range of 60 dB, quantised in Vi dB steps 
and with a resolution of 200 elements on the time 
axis. The display of the alphanumeric data for the 
annotation and the table of results was fitted into 
the remaining spaces. 

Apart from the memory, required to store 
the data for the 256 x 256 picture elements 
(8k bytes) and the program (2k bytes), the display 
unit also contains a single-chip display controller 
(Motorola MC 6845) and a microprocessor 
(Motorola MC 6802). This gives the display a 
computing ability which greatly simplifies the 
transfer of data to the display unit. The software 
for. the display includes a range of line, graph and 
block drawing facilities and a user-programmable 
character set. All of the interpolated points for the 
graphical functions are generated by the display 
software. 

The display format is a standard 625 iine/50 
field per second television raster scan which 
enables any standard picture monitor to be used 
directly. A 225mm diagonal picture monitor is 
built into the equipment and a separate video 
output socket is provided for driving an external 
monitor. Communication between the main 
controller microprocessor and the display unit is 
via a standard RS 232 serial interface which 
requires only two wires, thus simplifying the inter- 
connecting wiring. This interface also allows the 
display unit to be used as a graphics display for 
other purposes. 



5.6. Main controller 

The main controller is required to provide 
the sequence of instructions for the remaining 
units appropriate for the reverberation time (R.T.) 
range and frequency band settings. In addition, it 
has to interpret the inputs from the keyboard and 
the GPIB interface (see Appendices, Sections 12.3. 
and 12.4.). The total number of sequence 
variations possible is so large that a fixed 
programme sequencer would have been very 
complex. Therefore, a microprocessor is used to 
derive the appropriate sequence from a relatively 
small fixed programme memory. Once a micro- 
processor has been incorporated, a number of 
useful functions such as automatic computation of 
the slope of the decay and display of the 
reverberation time in numeric characters can be 
added at very low cost. 

The control system is based on an Intel 
8085A microprocessor. In addition to the system 
software which occupied almost 8k bytes of read- 
only-memory,- provision was made within the unit 
for 8k bytes of random access memory of which 
this equipment uses only 2k bytes. A 24-line 
parallel input/output device, a bi-directional serial 
interface, a counter/timer and a GPIB interface 
device 13 - 14 are also included in the design. The 
serial interface provides a low-speed system input 
for the keyboard and a high speed output for the 
display. The counter/timers are used to control 
the timing for the instrument when measuring 
events in real time (for example, the sound 
pressure level decay waveform). The GPIB 
interface, using the address set by switches 
connected to the parallel input/output device was 
provided to allow control of the instument and the 
processing of the results as part of a future, larger 
integrated acoustic measurement system and 
conformed to the IEC 625 and IEEE 488 standards 
for the interconnection of measuring instruments. 

The entire controller was designed on one 
printed circuit card and was intended as a small, 
self-contained microcomputer which could be 
used for other purposes besides the one described 
in this report (see Section 9). 

5.7. Keyboard 

The keyboard design is based on a National 
Semiconductor MM 57499 keyboard encoder which 
provides all of the waveforms necessary to scan the 
keyboard, interpret the inputs and send the 
appropriate encoded data via a serial output 
section to the main controller. A keyboard was 
built into the complete instrument and provision 
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also made for an additional keyboard which could 
be a small, hand-held remote control unit. 

5.8. Power Supplies 

The essential power requirements of the 
complete instrument totalled about 30W which 
made the option of battery operation feasible. 
However, because of the attendant requirement for 
an audio power amplifier with a maximum power 
ouput of about 50W, the combination would make 
battery operation for a worthwhile period of time 
impractical. This total dependency on a supply 
of mains power should not be a serious 
disadvantage in practice. 

The secondary power requirements of the 
instrument are 3 amps at 5 volts, 1 amp at 12 volts 
and 100 mA at both +15 and -15 volts. This 
corresponds to a total primary input power of 
about 46W which is low enough not to require 
forced ventilation of the instrument case. 

5.9. Prototype construction and enclosure* 

The printed circuit board designs are based 
on a BBC standard 4U high (144mm) card, 250mm 
long. Six boards were designed, one each for the 
microphone amplifier and band-pass filter, rectifier 
and analogue-to-digital converter, excitation 
generator, control system, display driver and the 
power supplies. A backplane was also designed to 
carry the sockets and to provide the 
interconnections between these six boards. These 
boards, together with the display monitor and a 
front panel carrying the keyboard are fitted into 
a standard 475mm wide x 4U high rack. 

The rear panel retains the boards in position 
and carries connectors for the input and output 
signals. Two audio frequency inputs to the micro- 
phone amplifier board are provided, one for 
microphone4evel signals and a second, which 
bypasses the microphone preamplifier, and which 
provides a high-level input. The excitation 
generator has one socket for the signal output to 
the loudspeaker amplifier and an input socket for 
an alternative excitation signal generator. Two 
multipin sockets are provided, one for the GPIB 
interface and printer and an auxiliary socket for 
other connections such as an external keyboard. 
A mains supply input socket and a mains power 
switch are also provided on the rear panel. 

The complete instrument is constructed 

*Most of the development of the prototype was carried out by 
K.J. Wright. 



in a proprietary frame which also serves as a 
carrying case. 

6. Software 

When the instrument is switched on the 
controller software first sets up the basic display, 
Fig. 6. The display area is split into four main 
parts. Two of these comprise the graphical display 
areas, each of which spans a 60 dB range of signal 




Fig. 6 - Initial display. 

levels on the vertical axis and contains 200 sample 
points along the horizontal axis. Both axes are 
subdivided by markers, the vertical axes into the 
equivalent of lOdB intervals and the horizontal 
axis into 5 equal spaces. The horizontal axis is 
also annotated to represent time, — 0.5 seconds 
in the example shown in Fig. 6. The third part 
of the display area, on the right hand side, consists 
of a table containing the ISO band numbers of the 
1/3-octave wide measurement frequency bands, 
starting at band 17,50 Hz, and the measured 
reverberation time in each band. The current 
operating frequency band is indicated by a 
highlight cursor. The fourth part of the display 
area consists of a single line of numbers and 
symbols at the top edge of the display area. These 
will be described below. The software also 
contains routines for recognising and interpreting 
the inputs from the control keyboard, an external 
keyboard and the GPIB interface. 

On receipt of the "start" instruction, the 
microprocessor sets the excitation generator and 
the microphone amplifier filter frequencies to the 
current operating frequency and the microphone 
amplifier gain to dB. Any traces remaining in the 
display areas from earlier measurements are erased. 
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The timing clock frequency is set to 200 times the 
reciprocal of the current range setting to give 200 
samples in the time period corresponding to the 
full scale range. The output from the excitation 
generator is then enabled for 115 clock periods*. 
During the first part (85 clock periods) of this 
time, the output of the analogue to digital 
converter is sampled once per clock period and the 
samples are integrated. Using these integrated 
samples, the additional microphone amplifier 
gain required to raise the mean level of the micro- 
phone signal to —5 dB relative to the full scale 
amplitude can be calculated. This gives an 
available margin to allow for signals which increase 
momentarily after the cessation of the loudspeaker 
output and also to accommodate a 5 dB ratio of 
peak-to-mean signal levels. The microphone 
amplifier gain is then increased by that amount 
and the excitation of the room continued for an 
additional 30 clock periods to allow the input 
signal processing circuits to settle at the new gain 
setting. 

At the end of the excitation period the 
output signal is terminated and the 'trace' mode is 
begun. To reduce the generation of unwanted 
modulation sidebands, the output of the excitation 
generator is not switched off abruptly but 
progressively reduced in amplitude at a rate of 
3 dB per step for the first 20 of the 200 clock 
periods of the 'trace' mode. From the beginning 
of the 'trace' mode, the output of the analogue 
to digital converter representing the microphone 
signal level is displayed as a function of time in 
the upper half of the graphical display area. The 
lower half of the graphical display area is used in 
a similar way to display the point-by-point average 
of several (or many) such traces. At the end of 
the first trace period, the trace number (T) in the 
top line of the display area is set to one. Fig. 7 
shows the appearance of the display after four 
traces. 

This cycle of excitation and trace continues 
automatically, the trace number being incremented 
by one at the end of each trace. The microphone 
amplifier gain set by the controller during each 
excitation period is displayed as the second 
number in the top line of the display (G). The 
averaging function, used to generate the average 
trace, is variable to give a range of exponential 
'time constants' from about 1.4 to 254 traces and 

*This time must be long enough to allow the second pressure level 
in the room being measured to reach its asymptotic value. The 1 15 
clock periods greatly exceed the time which is theoretically 
necessary but instrumental constraints make it difficuit to achieve a 
significantly shorter period. 




Fig. 7 - Display after several traces at one 
frequency. 

the multiplying factor is displayed as the third 
number (AC) in the top line of the display area. * 
The fourth number (N) is a counter which can be 
used automatically to terminate the sequence after 
a presettable number of traces. 

The main controller software also contains 
the computational algorithm** which uses the data 
stored for the average decay trace to compute the 
mean slope of the decay and, together with the 
current range setting, the corresponding 
reverberation time. It also calculates the position 
of the straight line corresponding to the 
reverberation time and superimposes this line on 
the display of the average trace. At the same time, 
the entry in the table of reverberation times at the 
current operating frequency is changed to the 
new calculated value in seconds, as depicted in 
Fig. 8. This routine can be initiated manually 
at any time by the operation of the appropriate 
front-panel key or automatically when the total 
number of traces executed is equal to the number 
set in the preset counter. In some circumstances, 
the computed straight line is not the best fit 
to the data and provision has been made for 
incrementing or decrementing the result under 
manual control. These circumstances can arise if 

*The average trace is generated by multiplying the data for the 
previous average trace by the averaging constant, AC, (which is in 
the range 0.5 —1.0), adding the new trace data and then scaling the 
sum to normalise the result. This process implements an 
exponential decay in the domain of the traces. The equivalent 
exponential time constant is equal to — 1 /in (AC) traces. 
**This algorithm is described in more detail in Section 7. 
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Fig. 8 - Display after computation of reverberation 

time. 

the average decay is significantly curved or if the 
dynamic range of the measurement is inadequate. 
A human operator can often make better 
judgements in these cases. It is probable that the 
algorithm could be improved to give more accurate 
results for these difficult cases. 

Once the measurement is completed at one 
frequency, the frequency cursor can be moved in 
either direction to repeat the measurements at 
other frequencies. Provision was made in the 
software for completely unattended operation; 
after executing the preset number of traces at one 
frequency, the controller automatically steps on to 
the next highest frequency and repeats the 
measurement procedure until it has completed the 
measurement in the highest frequency band. Both 
the preset number of traces and the time ranges 
are independently presettable at each frequency. 

The controller software also contains 
routines for remote operation by an external 
master controller via the GPIB Interface. This 
interface system allows the interconnection of up 
to 15 suitably equipped separate instruments to 
form a larger measurement system under the 
control of the master controller. All of the 
keyboard functions can be implemented via the 
GPIB and, in addition, the computed results and 
the data representing the individual and the average 
traces can be transmitted via the GPIB back to the 
master controller. Provision in the software is also 
made for the connection of a simple, compatible 
printer to the GPIB terminals to print out the 



reverberation-time table directly, without the use 
of a master controller, using the 'TALK ONLY' 
mode of the instrument's GPIB interface. 

Fig. 9 shows a drawing of the prototype 
keyboard layout and the Appendix (Section 12.4) 
contains an annotated list of the keyboard and 
GPIB instrument functions. 

7. Computation algorithm 

An idealised trace of the sound pressure 
level in decibels as a function of time, beginning 
at the time the excitation of the room ceases, 
consists of three straight line portions. The first 
portion, of zero gradient, represents the time 
delay of the sound energy from loudspeaker 
to microphone. The third portion, also of zero 
gradient, represents the constant background 
noise level in the room. Between these two 
portions is the line of finite gradient representing 
the decay of the sound energy as a function of 
time. The gradient of this middle portion, when 
appropriately scaled, is equal to the reciprocal of 
the reverberation time. In practice, real traces 
differ from the ideal in having deviations from the 
ideal trace which are sometimes very large. Some 
of these deviations are variable in that they differ 
from trace to trace and are therefore reduced in 
effect by the trace averaging process. The other 
deviations are fixed from trace to trace and are 
therefore not reduced by the trace averaging 
process. These latter deviations indicate the lack 
of perfect diffusion of the sound energy within the 
room 15 or the presence of multiple resonating 
acoustic or mechanical systems. 

The main problem in extracting the gradient 
of the middle portion of the composite trace is the 
location of the two transition points between the 
three portions of the trace, especially in the presence 
of irregular deviations from the ideal trace. Because 
these transitions are not abrupt, even for an ideal 
trace, the calculation of the gradient of the line has 
to be confined to that part of the middle portion of 
the trace which is not within about 5 dB of either of 
the transition points. Five decibels represents an 
approximate threshold beyond which one portion of 
the line has very little effect on the adjacent portion. 

The gradient of the middle portion of the 
line can be satisfactorily estimated using linear 
regression by calculating the position and slope of 
the straight line which has the least total mean 
squared difference between itself and the data 
points representing the actual sound energy decay. 
A weighting function was introduced into the 
mean squared error calculation which weighted the 
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Fig. 9 - Prototype keyboard layout. 



data according to the local average slope of the 
trace, local being ± 5% of the whole trace duration. 
In this way, by choosing a suitable weighting 
function, the data representing the first and last 
portions of the trace, having theoretically a slope 
of zero, are assigned zero or a very small weight 
whereas the middle part of the trace, having 
a finite slope throughout, is assigned a high 
weighting and thus contributes most to the mean 
squares calculation. 

In practice, this algorithm works very well 
even with moderately irregular traces. However, 
extensive comparison with other methods of 
measurement showed that the computed 
reverberation times were approximately 2% too 
high. This error was attributed to the fact that the 
third portion of the trace, that representing the 
background noise level in the room being 
measured, had, even after averaging several traces, 
sufficient irregularities to make the weighting 
coefficients non-zero. The gradient of the 
computed straight line was thus being reduced by a 
disproportionate amount by these components 
near to the end of the trace. An additional 
function was therefore introduced into the 
algorithm to eliminate entirely that portion of the 
trace which was close to the background noise 
level. A sample at the end of the trace is taken to 
represent the background noise level and all data 
which does not exceed this threshold by at least 
6 dB are assigned a weight of zero. An Appendix 
(Section 12.5) gives a full description of the final 
version of this algorithm. 

Tests on the modified algorithm showed no 
measurable consistent error for traces which were 



reasonably free of regular departures from a linear 
decay. In this form, the algorithm is still incapable 
of correctly interpreting traces with pronounced 
curvature, and also gives inaccurate results when the 
signal-to-noise ratio is less than about 10 dB. Despite 
the large amount of arithmetic involved in 
calculating the local slope over 19 sample points 
and the sums and products required for the mean 
squares fitting at every one of 180 sample points, 
the computation and display of the result requires 
only about 0.4 seconds, which is sufficiently fast 
not to constitute a serious delay in the room 
measurement procedure. 

The algorithm would have to be improved to 
give an automatic indication of failure for 
significantly curved or irregular decays before 
the instrument would be capable of completely 
unattended operation in any area. However, 
the requirement for unattended operation with this 
kind of instrument only applies for the 
measurement of materials in a reverberation 
room where such curved decays should not be 
present. Alternatively, most of the curved decays 
observed in studios and control rooms are present 
only at one of the usual five microphone positions 
and could be reduced significantly by including in 
the trace averaging process a space-average by using 
either several microphone positions* or a rotating 
boom microphone. Work continues on this aspect 
of reverberation time measurement, but no 
significant results were available at the time of 
writing. 

*Several microphone positions are always used in reverberation 
time measurements but the results are normally noted individually 
and averaged arithmetically to obtain the space average. 
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8. Operational tests 

After the completion of the prototype 
equipment, Fig. 10, it was used on -..several 
occasions for field measurements. Although it was 
fairly delicate and suffered from a number of 
minor faults, most of these faults were found to be 
essentially mechanical and due to the prototype 
construction. No significant faults were found 




Fig. 10 - Completed prototype instrument. 

in the design or in the operational principles 
although the initial keyboard layout was incon- 
venient and was improved for the final de'sign. 
Changes were made during the course of the 
operational tests to the gain range and frequency 
response of the microphone amplifier, to reduce 
the incidence of amplifier overload by out-of-band 
signals. The gain range had originally been 52 to 
112 dB for all frequencies from 10 Hz to 15 kHz. 
It was possible for the microphone amplifier to 
overload on the low-frequency background noise 
level when attempting to measure high-frequency 
reverberation times with a low peak sound pressure 
level (i.e. when trying to avoid disturbing areas 
adjacent to the one being measured). The 
amplifier gain range was reduced to 40 to 100 dB, 
which was still sufficient to show the inherent 
microphone amplifier noise on the display. A 
single-pole, high-pass filter giving 18 dB of 
attenuation at 50 Hz was also inserted between the 
microphone pre- and main amplifiers. 

Fig. 11 shows the appearance of the 
display after a complete measurement at one 
microphone position. The reverberation time 
measured at each frequency has been entered 
and stored in the tabulated results ready for 
noting, printing, or transmission to the 
controlling computer. 




Fig. 11 - Display after complete measurement at 
one microphone position. 

9. Future applications 

The complete instrument for measuring 
reverberation time may be used for other types of 
acoustic measurements with changes only to the 
control software. These include the measurement 
of sound insulation, loudspeaker power output and 
sound-pressure-ievel data logging. 

There is, however, one limitation to the use 
for these other purposes and that is the accuracy of 
the digitally controlled attenuator which sets the 
microphone amplifier gain. The device used had a 
specified accuracy of ±3/4 dB which is not 
sufficiently accurate for precision laboratory 
acoustic measurements, except reverberation time 
for which the absolute gain is unimportant. This 
problem could be overcome if a more accurate 
device were available, or if the measurement was 
such that the gain need not be changed. 
Alternatively, a separate initial calibration could be 
carried out and the corrections permanently stored 
as functions of frequency and amplitude. 

The individual units from which the instru- 
ment is made may also be used for other purposes. 
In particular, the display controller was specifically 
designed to be applied generally as an intelligent 
display with both alpha-numeric and graphics 
facilities. The main controller was also specifically 
designed as a small, single-board, general-purpose 
microcomputer and would be capable of being 
used in a wide range of applications. It appears to 
be unique amongst small, single-board micro- 
computers in that it includes a complete GPIB 
interface. 
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10. Conclusions 

The design of prototype new equipment for 
measuring reverberation times in a wide range of 
different rooms has been successfully completed. 
The new equipment includes all of the functions 
of the equipment which it is to replace plus some 
additional functions which will simplify the 
routine measurement of both room acoustics and 
materials. It is also much lighter and requires only 
one person to operate it effectively instead of the 
two required previously. The final component cost 
(« £700 at 1981 prices) has confirmed the initial 
assertion that additional copies could be built at a 
low unit cost. 

The provisions made for remote operation 
will allow the instrument to be used as part of a 
future, integrated acoustics measurements system 
which would have completely automated 
measurement and data logging facilities. However, 
some improvements would have to made to the 
reverberation-time computation algorithm before 
completely unattended operation would be feasible. 

The individual component units of the 
instrument may be applied in other ways; the 
designs make their inclusion into other micro- 
processor systems a simple matter. 
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12. Appendices 
12.1. A cost-benefit analysis for the choice of the location of the analogue-to-digitai convertor ' a ' 



a.d.c. location 


Consequences on ... 


a.d.c. 


microphone amplifier 


filter 


logarithm & 
spl converter 


1 . After 

microphone 


impractical 
because of 
low signal 
levels 


not required 


complex, 

hi S h Lb) 

speed 
digital 


simple 
digital 
averaging, 
logarithmic 
conversion by 
look-up table 


2. After 

microphone 
amplifier 
(not gain 
controlled) 


very wide 

dynamic 

range 

> 16 bits, (b) 

highspeed 


analogue, simple 
but wide dynamic 
range 


3 . After gain 
controlled 
microphone 
amplifier 


simple, <12 
bits, high 
speed (b ' 


analogue, gain 
controlled 


4. After 
filter 


Analogue 

with 

switchable 

centre 

frequency 


Analogue 
rectifier and 
integrator. 


5 . After 
logarithmic 
and spl 
converter 


very simple, 
8 bit, low 
speed (c) 


Rectifier 

moderately 

difficult (d) 



Notes: (a) Analysis based on 60 dB display range of sound pressure level, 100 dB input dynamic 
range and 24 x l/3rd octave frequency bands (50 Hz to 10 kHz), (b) High speed denotes sampling the 
input audio signal at greater than twice the maximum signal frequency, that is > 25 kHz sampling 
frequency, (c) Low speed denotes sampling of the sound pressure level as a function of time (< 400 Hz 
sampling frequency), (d) In the final design, the logarithmic converter was included as part of the analogue- 
to-digital converter at a negligible additional cost. 

12.2. Giossary of terms 

Some of the terms used which are not common in acoustics are defined in the context of this 
report. 

active filter — a frequency selective network containing amplifying or impedance-changing components. 

address — a number specifying the location within a memory space of an element of data. 
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alpha-numeric 

analogue 
binary divider 

byte 
Chebyshev filter 

elliptic filter 
field scan 

GPIB 

input offset 
IEC 625 

IEEE 488 

Kbyte or K 

line scan 

maximal length 
sequence generator 

microprocessor 

peripheral 
phantom power 



random-access 
memory (RAM) 

raster scan 



read-only 
memory (ROM) 



printed or displayed data using the alphabetical letters, numerals and 
punctuation marks only 

representation of a quantity or value by a continuous function 

a system for reducing the frequency of an alternating signal by successive 
halving 

a binary coded word of 8 digits 

a class of filters, based on Chebyshev polynomials, in which the zeroes of 
the transfer function are distributed throughout the passband in such a 
way as to equalise the passband ripple 

a class of filters, based on elliptic functions, in which the poles and zeroes 
of the transfer function are mutually reciprocal — first introduced by Cauer 

the lower of the two scanning frequencies of a raster scan 
display 

General Purpose Interface Bus (see IEC 625 and IEEE 488) 

in differential input amplifiers, the input required to produce zero output 

International Electrotechnical Commission standard for the inter- 
connection of measuring equipment 

Institute of Electrical and Electronic Engineers standard for the inter- 
connection of measuring equipment 

common abbreviation for 2 10 (1024) bytes 

the higher of the two scanning frequencies of a raster scan display (q.v.) 

an N-stage register (q.v.) which has a feedback connection that results in 
an output sequence of length 2 N —1 periods 

the essential elements of computer central processor integrated on a single 
chip 

external to the main part of the computer 

operating power supplied to a microphone amplifier via the signal 
frequency connections 

a computer data store in which the data may be accessed in random 
order for writing or reading variable data 

a method of data display in which the display area is periodally described 
by two orthogonal deflection systems 

a computer data store in which the data is fixed by a process outside 
the computer system — usually used for fixed programs 
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RS 232 — EIA standard for implementation of a serial interface 

shift register — a connection of data storage elements in which the data passes 

sequentially from one location to the following location 

TALK ONLY — a mode of the GPIB interface in which the instrument autonomously 

transmits data without the intervention of a master controller 

12.3. General purpose interface bus (GP!B) 

12.3.1. General 

This bus interconnection system is based on developments by the Hewlett-Packard Company which 
have been established through the IEEE as an American industry standard (IEEE 488) and by the IEC 
as an international standard (IEC 625) for interconnecting up to 15 digital instruments in order that they 
may exchange data. Each of these instruments must of course conform to the standard. The system is 
fairly complicated and no more than a brief outline can be given here. Names and expressions which have 
a specific function, defined in the GPIB standards are given here in the upper case. 

Each of the subordinate instruments may either send (TALK) or receive (LISTEN) as appropriate 
and the system must contain at least one master instrument (CONTROLLER), which is usually either a 
computer or a programmable calculator, to organise the system operation. Only one TALKER is permitted 
at a time but any number of instruments may LISTEN simultaneously. Commands issued by the 
CONTROLLER assign individual instruments to be either TALKERS or LISTENERS. 

In very small systems, individual instruments may be able, if suitably programmed, to either TALK- 
ONLY or LISTEN-ONLY and under these circumstances no CONTROLLER is required. For example, a 
single instrument (set to TALK ONLY) may send data to a printer (set to LISTEN ONLY). 

The interconnection system consists of a 24 wire bus which can be functionally divided into 3 groups. 
Eight wires, DI01 - DI08 make up the data group which transmits the data (in 8-bit bytes) in either 
direction. A further eight wires make up the command group and are used to synchronise the transfer of 
data, to differentiate between data and instructions on the data group and also to transmit general 
commands. The remaining eight wires are signal returns. 

The functions of the 16 active lines are defined as follows: 

Data Group 

DI01 — DI08 — Data inputs/outputs 1 to 8, used for bi-directional interchange of data 

and multi-bit commands 

Command Group 

DAV — Data valid, issued by the TALKER to indicate that new data is available 

on DI01 - DI08 

NRFD — Not ready for data, issued by LISTENERS to indicate that they are not 

yet ready to receive new data 

NDAC — Data not accepted, issued by LISTENERS to indicate that they have not 

yet finished with the data currently on DI01 - DI08 

IFC — Interface clear, issued by the CONTROLLER to set all instruments to an 

idle state 
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ATN - Attention, issued by the CONTROLLER to identify commands on DI01 

-DI08 

SRQ — Service request, issued by an instrument to indicated that it has completed 

a task and/or requires further instructions 

REN — Remote enable — issued by the CONTROLLER to enable instruments to 

switch to remote control 

EOI — End or identify, issued by a TALKER to indicate the end of transmission 

sequence or to identify the instrument which sent SRQ. 

The first three lines of the command group, DAV, NRFD and NDAC, control the timing of the data 
exchange in such a way that the rate of data transfer is set by the slowest instrument actively connected to 
the bus. For example, each LISTENER holds NDAC true until it has accepted the data or command 
currently on DI01 — 8. Only when all LISTENERS have accepted the information does NDAC go to 
false, allowing the TALKER to remove or replace the information. A similar process occurs for NRFD. 

On completion of a task, an instrument which might not be enabled either as a TALKER or as a 
LISTENER may issue a "request for service" on the SRQ line. The subsequent sequence of events depends 
on the programme of the CONTROLLER but, usually, the controller begins either a SERIAL POLL or a 
PARALLEL POLL in order to find out which instrument issued the request. In a SERIAL POLL each 
instrument is enabled as a TALKER in turn and a status word and an output on the EOI line is returned by 
the instrument which made the request. In the PARALLEL POLL mode, up to eight instruments are 
enabled simultaneously and each instrument replies using one preselected line of the data group. 

Apart from the special case of the PARALLEL POLL mode, instruments are normally identified by a 
5 -digit binary code word which is preset in the instruments, normally by internal switches or pre-wired 
connections. By sending a command consisting of the appropriate five digit code, together with a further 
two digits to specify the mode, any individual instrument may be set to one of four different modes. 
Three of these modes are TALK, LISTEN and SECONDARY COMMANDS. One of the five digit binary 
code words is reserved to implement the general UNTALK and UNLISTEN commands. The fourth of the 
available modes is used to identify other general commands. In addition, if any instrument which is in the 
TALK state receives a TALK address which is not its own then it will also go to the UNTALK state. 

12.3.2. Implementation in the reverberation time measuring equipment 

The majority of functions defined in the standards are implemented in the reverberation time 
measuring equipment. Those not included are EXTENDED ADDRESSING, PARALLEL POLL, DEVICE 
CLEAR and DEVICE TRIGGER. 

The function subsets (as defined in reference 13) are 

SHI — source handshake (talker) 

AH1 — acceptor handshake (listener) 

T5 — talker (full capability) 

L4 — listener (but not listen only) 

SRI — service request 

RL1 — remote control 

PPO — parallel poll, not implemented 

DCO — device clear, not implemented 
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DTO — device trigger, not implemented 

CO — controller functions, not implemented 

12.4. Instrument functions 

The functions available in the prototype reverberation time measuring equipment are initiated by 
transmitting a single ASCII coded character to the instrument, either from the internal keyboard or via the 
GPIB. Although the internal keyboard actually sends single characters, the keys are labelled with the 
function name. The functions incorporated (together with the key label) are 

A — Increment current frequency (CURSOR | ) 

B — Decrement current frequency (CURSOR^) 

C — Increment multiplying constant for trace averaging (AVERAGING TIME 

CONSTANT 4 ) 

D - Decrement multiplying constant (AVERGING TIME CONSTANT | ) 

E — Increment time scale (RANGE 4 ) 

F — Decrement time scale (RANGE f ) 

G - Increment preset counter (NO. OF TRACES 4 ) 

H — Decrement preset counter (NO. OF TRACES ? ) 

I — Start (erases previous trace data) (START) 

J - Stop (STOP) 

K — Continue (does not erase previous trace data) (CONTINUE) 

L - Compute and display reverberation time (REVERBERATION TIME 

COMPUTE) 

M — Start automatic mode. Instrument executes traces until preset count is 

reached, computes reverberation time then steps to next highest 
frequency and continues measurement (AUTO) 

N — Start semi-automatic mode. Same as automatic mode but instrument 

stops instead of stepping to next highest frequency (SEMI AUTO) 

-*■ - Increment calculated reverberation time (REVERBERATION TIMEf) 

-*- - Decrement calculated reverberation time (REVERBERATION TIME f ) 

X — Shift key which can be operated once before instructions E, F, G and H 

to make these four functions affect all 24 frequency bands instead of just 
the current operating frequency. This shift function is cancelled by 
operation of any key other than E, F, G and H. An indicator light is 
included to show when this function is active (ALL BANDS) 

In addition to the above, several extra functions are included which are intended to be initiated via 
the GPIB. These are -.- 

O — Current operating frequency band set by the next two decimal digits 
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Q 



R 



U 



transmitted in the range 17—40 

— The averaging constant is set by the next character transmitted according 
to Table 1. 

— The time scale is set by the next character transmitted according to 
Table 1 

— The preset counter is set by the next two decimal digits transmitted in 
the range 00 — 99 

— The current measured reverberation time is returned via the GPIB to 
any listeners 

— The 200 data samples representing the current trace are returned in four 
blocks of 50 samples* via the GPIB to any listeners 

— The 200 data samples representing the average trace are returned in four 
blocks of 50 samples* via the GPIB to any listeners 

Table 1 - Second Character Function for Functions P and Q. 



Second 
Character 


Function 


Set AC 
(First character-P) 


Set RANGE 
(First character-Q) 


A 
B 
C 
D 
E 
F 
G 
H 
I 

J 
K 













1 


5 

75 

80 

85 

9 

92 

94 

96 

98 
99 

00 


0-.5 s 
0-.75 s 
0-1.0 s 
0-1.5 s 
0-2.0 s 
0-2.5 s 
0-3.75 s 
0-5.0 s 
0-7.5 s 
0-10.0 s 
0- 15.0 s 



12.5. Reverberation time computation algorithm 

The data on which the computation of the reverberation time is based are those representing the 
averaged sound pressured level decay as a function of time. The data is stored as 200 eight-bit binary- 
coded samples, each representing an averaged sound pressure level in the range to —63.75 dB (the negative 
sign is not stored but assumed — it would otherwise have taken an additional bit). The signal representation 
is in the form of sign and magnitude, —63.75 dB being represented by the binary-coded equivalent of 
decimal 255. With eight bits available to represent the whole signal range, each quantisation level represents 
KdB. 

The algorithm can be divided into four distinct stages, the assessment of the background noise level, 
the computation of the weighting function, the summation of the partial products for the mean square 
error calculation and, finally, the calculation of the slope and the intercepts with the axes of the straight 
line representing the best fit with the real data. 

*The data array has to be subdivided because of the limits on block length inherent in most GPIB controllers. These vary between 64 and 80 
characters, depending on the controller. 



(PH-235) 



- 21 



12.5.1. Assessment of the background noise level 

The background noise level in the room being measured or the level of the decaying sound at the end 
of the measurement time period if the decay is not complete by that time, is simply assessed by taking one 
sample at the end of the trace (sample number 199) and storing this value minus 24 (+6 dB) as the 

threshold. 

12.5.2. Computation of the weighting function 

The weighting function required for the elimination of the first and third parts of the trace from the 
computation is calculated from the local slope at each sample position, beginning with sample number 10 
and ending with sample number 189. For each sample, the sum of the 9 preceding sample values is 
subtracted from the sum of the 9 following sample values to give a number which is positive for sound 
pressure levels which decrease with increasing time. 

Including the negative sign preceding the level in dB, for a uniform gradient of — s dB/sample 
(= — 4 s quantisation levels/sample) and a middle sample value of —X dB, this calculation can be 

represented by:- 9 9 

F = 2 (-4X + 4xh)-2 (-4X - 4 s n) 
n=l n=l 

where F is the positive numerical value representing the local gradient in terms of quantization 
levels/sample 



This can be simplified to 



or 



9 

F = 8 s 2 (n) 

1 

Y=360s 



The weighting coefficient, w, used for each sample is assigned according to the following table, where 
approximate decimal equivalents to the binary-coded fractions are given. 

w = for F<50 

w = 0.004 for 50<F<75 

w = 0.008 for 75<F<100 

w = 0.016 for 100<F<125 

w = 0.031 for 125<F<150 

to = 0.063 for 150<F<175 

w = 0.125 for 175 < F < 200 

w = 0.250 for 200<F<225 

w = 0.500 for 225 < F < 250 

w = 0.996 for 250 < F 

The weighting coefficient increases in value rapidly for values of F greater than about 200, corres- 
ponding to approximately 0.56 dB/sample or the complete 60 dB of displayed signal range in 108 sample 
intervals. The effect of this weighting coefficient for different gradients can be illustrated graphically by 
plotting the weighting coefficient as a function of the number of sample periods which the local gradient 
would require to cover the whole of the 60 dB displayed range if that gradient was continued over the 
whole range. Fig. 12 shows such a graph for the prototype weighting function. It is probable that this 
function could be improved; for example, a regular decay which required 200 sample periods to cover the 
60 dB range would be assigned a low weighting coefficient for most sample values and the result would 
therefore be unduly affected by any local departures from the average line. 

12.5.3. Computation of mean square error 

This computation is carried out for each of the 180 samples in turn, immediately following the 

(PH-235) — 22 — 



1-0 



g0-8 



§ 0-6 
o 



<?0-4 



fo-z 



limit of 
displayed 
samples 

I 

I 

I 

I 

I 



iOO 200 300 

number of sample periods 
required to cover 60dB range 



Fig. 12 - Weighting coefficient function. 



computation of the weighting coefficient described in section 12.5.2. The five different summations 
required by the least mean square error computation are progressively accumulated. 

The first step in the computation of these summations is the check that the current sample value is 
less ,than (higher sound pressure level) the threshold representing the background noise level + 6 dB, 
obtained in stage 1 (12.5.1.). If the sample sound pressure level is lower than the threshold then the 
weighting coefficient, w, is set to zero and the sample thus excluded from the summations. 

The gradient of the best-fit straight line, g, is calculated from 16 * 

2 (w . x) . 2 (w . y) — 2 (w) . 2 (w . x . y) 



(2 (to . x) ) 2 - 2 (w . x 2 ) . 2 (w) 
and the intercept with the vertical axis, i from 

2 (to . x) . 2 (w . x . y) — 2 (w . x 2 ) . 2 (w . y) 



(X(w . x)) 2 - X(w . x 2 ) . X(w) 
where x represents the sample number 

y represents the sample value 
and all of the summations are carried out over the range of x from 10 to 189. 

12.5.4. Computation of the dispiayed results 

The result of the calculation of section 12.5.3 . is the equation of the straight line which gives minimum 
total square error between itself and the weighted original data. From this equation, the position of the 
line on the display screen can be calculated and the line displayed. From the gradient and the current range 
setting, the numerical value of the reverberation time can also be calculated and displayed. 

Rate of decay of sound energy, R = g quantisation levels per sample period. 



*The expressions for the slope and offset are the conventional expressions for the linear regression of y on X with the addition of the weighting 
coefficients, W. 
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Because there are four quantisation levels per decibel (Section 12.5.), 

R = g/4 dB/sample period. 
Also, as there are 200 sample periods in the time corresponding to the range setting, that is:- 

1 sample period = range/200 seconds (Section 5.5.) 
then 

R = 200g /(4. range) dB/s 

And the reverberation time, RT, that is the time taken for the sound energy to decay by 60 dB, is 
given by:- 

RT = 60 /R 

RT = 1.2 range /g 
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